[C n] Emission From NGC 4038/39 (The "Antennae") 
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ABSTRACT 



We present observations of NGC 4038/39 in the [C II] 158 /xm fine struc- 
ture line taken with the MPE/UCB Far-infrared Imaging Fabry-Perot Inter- 
ferometer (FIFI) on the KAO. A fully sampled map of the galaxy pair (with- 
out the tidal tails) at 55" resolution has been obtained. The [C II] emission 
line is detected from the entire galaxy pair and peaks at the interaction zone. 
The total [C II] luminosity of the Antennae is Lrcri] = 3.7 x 10 8 Lq, which is 
about 1% of the far-infrared luminosity observed with IRAS. The main part 
of the [C II] emission probably arises from photodissociation regions (PDRs), 
and a minor fraction may be emitted from H II regions. A small part of 
the [C II] emission comes from standard cold neutral medium (CNM); how- 
ever, for high temperature (T ~ 100 K) and high density (tih ~ 200 cm -3 ) 
about one third of the observed [C II] emission may originate from CNM. 
From PDR models we derive densities of the order of ~ 10 5 cm -3 and far-UV 
(FUV) intensities of 460xo, 500^0, and 240xo for the PDRs in the interaction 
zone, NGC 4038, and NGC 4039, respectively. However, PDRs with densities 
of the order of ~ 10 2 cm -3 and FUV intensities of the order of ~ 100xo could 
also explain the observed [C II] emission. The minimum masses in the [C II] 
emitting regions in the interaction zone and the nuclei are a few x 10 7 M Q . A 
comparison with single dish CO observations of the Antennae shows a [C II] 
to CO intensity ratio at the interaction zone a factor of 2.6 lower than usually 
observed in starburst galaxies, but still a factor of about 1.3 to 1.4 higher 
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than that at the nuclei of NGC 4038/39. Therefore, no global starburst is 
taking place in the Antennae. [C II] emission arising partly from confined 
starburst regions and partly from surrounding quiescent clouds could explain 
the observed [C II] radiation at the interaction zone and the nuclei, though 
the star formation activity toward the nuclei is lower. Accordingly there are 
small confined regions with high star formation activity in the interaction 
zone and with a lower star formation activity in the nuclei. This supports 
the high density and high FUV intensity for the PDRs in the interaction zone 
and the nuclei. 

Subject headings: galaxies: interaction — galaxies: individual: (NGC 4038/39) 
— infrared: galaxies — infrared: spectra 

1. Introduction 

The galaxy pair NGC 4038/39 (Arp 244) is an interacting system in an early stage 
of merging at a distance of about 21 Mpc from our own galaxy. On long-exposed im- 
ages in the optical (e.g. Arp 1966) the interaction is clearly visible because of the tails 
( "Antennae" ) emerging from two uniformly luminous, partly overlapping ovals and be- 
cause of the dwarf galaxy that appears to have formed at the tip of the southern tail 
through the interaction (Zwicky 1956, Schweizer 1978, Mirabel, Dottori, & Lutz 1992). 
The tails contain about 70 % of the total amount of H I in the system (van der Hulst 
1979). Short-exposed images (Laustsen, Madsen, & West 1987) reveal hints of the in- 
teraction on somewhat smaller scales; the distorted arrangement of Ha knots and the 
velocity distribution of the individual knots lead Rubin et al. (1970) to conclude there 
is an interaction of two rotating galaxies. Computer simulations carried out in the clas- 
sical paper of Toomre & Toomre (1972) and later by Barnes (1988) can account for the 
present morphological appearance of the system quite well by assuming an interaction 
of two rotating spiral galaxies. 

Spectra of the nuclei of both galaxies taken in the optical range do not look like pure 
starburst spectra but consist of a composition of early-type stars and late giants (Keel 
et al. 1985). The detection of bright near-infrared peaks at the nuclei lead Bushouse & 
Werner (1990) to the same result. Their images of the Antennae in J- and R-band and in 
Ha show also the same pattern of bright knots in the surroundings of NGC 4038 and in 
the bridge connecting both galaxies. The Antennae system as a whole shows a relatively 
low star forming efficiency according to the ratio Lir/M(H 2 ) ~ 8.55 measured by Young 
et al. (1986); it is only a factor of 3 higher than in the Milky Way. The ratio determined 
by Sanders & Mirabel (1985) is almost twice as high. However, they observed a smaller 
region in CO and therefore probably underestimated the molecular mass. Comparison of 
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the Antennae with the sample of interacting and isolated galaxies of Young et al. (1986) 
shows NGC 4038/39 to have characteristics more like isolated galaxies. 

Measurements in the radio continuum at 1.5 GHz and 4.9 GHz of NGC 4038/39 
(Hummel & van der Hulst 1986) reveal a number of discrete knots which coincide in 
general with Ha knots, and an underlying diffuse component. This diffuse component 
has a steep spectral index on average which indicates non-thermal emission, and the 
peak of the diffuse radio emission is at the dust patch near the overlapping region. The 
discrete radio knots account for roughly 35 % of the total radio emission and have a 
spectral index of a ~ —0.5 on average, probably due to a thermal contribution (Hummel 
& van der Hulst 1986). 

Interferometric observations of the Antennae in CO (1 — > 0) by Stanford et al. 
(1990) show three main concentrations of CO emission. Two are associated with the 
nuclei and the third with the interaction zone. The overlap region is the strongest CO 
source and contains xs 10 9 Mq of gas, roughly as much H 2 as both nuclei together. Based 
on 10/im and Ha data, the authors have calculated a star forming rate of 5Mq yr _1 , and 
consequently the life time of the molecular gas is 2 x 10 8 yr. Single dish observations in 
CO were made at the interaction zone and both nuclei of the Antennae by Aalto et al. 
(1995). The ratio of the emission lines of 12 CO and 13 CO measured in the nuclei and 
the overlapping region of the Antennae is similar to that found in the central regions of 
"normal" starburst galaxies (Aalto et al. 1995). 

An excellent tracer of star formation activity in galaxies is the strong [C II] 158 jum 
2 P3/2 2 Pi/2 fine structure line which arises mainly from photodissociation regions 
(PDRs) created by far-ultraviolet photons from hot young stars impinging on nearby 
dense interstellar clouds (Crawford et al. 1985, Stacey et al. 1991). In combination with 
CO and FIR observations, the [C II] emission can be used with PDR models (Tielens 
k Hollenbach 1985, Wolfire, Hollenbach, & Tielens 1989, Wolfire, Tielens, & Hollenbach 
1990) to derive densities and far-UV intensities and estimates of the star formation 
activity. Extragalactic surveys of [C II] emission (Crawford et al. 1985, Stacey et a/.1991) 
concentrated mainly on nuclei while more recent observations have imaged individual 
galaxies to study the distribution of [C II] emission on large scales. [C II] images of M83 
(Geis et al. 1998) and NGC 6946 (Madden et al. 1993) demonstrate that the emission is 
extended at least over the full optical extent and often follows the distribution of the FIR 
and CO within the disk of the galaxies. In the case of NGC 6946 [C II] emission beyond 
the optical extent of the galaxy has been found. This [C II] emission has been attributed 
to diffuse gas and not to PDRs. Therefore we also estimate the possible contribution of 
[C II] emission from neutral atomic gas and from ionized gas in NGC 4038/39. Because of 
its relative proximity, NGC 4038/39 is a unique source for carrying out spatially resolved 
measurements of the [C II] line in an interacting system. We present the results of our 
imaging spectroscopy study of the [C II] line in NGC 4038/39 and compare them with 
observations obtained with ISO. 
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2. Observations 

The observations of the Antennae in the [C II] 2 P 3 / 2 — > 2 Pi/2 fine structure line at 
157.7409 /im have been carried out during four individual flights in 1992 with the Kuiper 
Airborne Observatory (KAO) from Christchurch, New Zealand, using the MPE/UCB 
Far-infrared Imaging Fabry- Perot Interferometer (FIFI; Poglitsch et al. 1991, Stacey 
et al. 1992). At 158 ^m the FWHM of the beam is about 55". The detector array 
consists of 5 x 5 pixels of size 40" x 40" per pixel on the sky. We observed three array 
positions to cover an area of 300" x 300" including both galaxies and the interaction 
zone and to fully sample an area of 150" x 150" around the interaction zone. Due to 
the large velocity dispersion in one array setting, as determined from the emission of 
the Ha knots (Rubin, Ford, & D'Odorico 1970) and from the CO emission (Stanford 
et al. 1990), we chose a spectral resolution of 144 km s _1 (FWHM) and a scan width 
of 380 km s -1 . The scan center was set at 1600 km s -1 . The wavelength calibration 
was implemented by means of the H 2 S absorption line at 157.7726 /im Flat fielding 
of the detector was carried out with two internal black bodies. The secondary mirror 
of the telescope was chopped approximately 4' in roughly the east-west direction, and 
the telescope was nodded to compensate for the beam offset. For absolute intensity 
calibration we observed Jupiter at 158 /im and assumed a temperature of 128 K at 
this wavelength (Hildebrand et al. 1985), an equatorial diameter of 43.22", and a pole 
diameter of 40.42" (The Astronomical Almanac 1992) for Jupiter. The accuracy of the 
intensity calibration is estimated to be about 30 % and the absolute pointing positions 
are uncertain by about 15". 



3. Data Reduction 

The Fabry- Perot interferometers in FIFI are adjusted to the appropriate velocity 
range and spectral resolution for every object. Further the Fabry-Perot interferometers 
are put to a 'park position' when the instrument is not in use. Therefore the Fabry-Perot 
interferometers are adjusted newly for every observing flight and for every object and 
there may occure slight shifts (~ a few km/s) of the velocity center of the spectral scans 
between individual observations. To combine the different pointings of the observation of 
the Antennae obtained in different observing flights we used the velocity range contained 
in all observations (1420 - 1740 km s _1 ). 

To combine the different array positions, to account for the statistical errors in the 
data, and to regrid the data points on a regular grid, we constructed a program using 
the method of maximum entropy. The program follows mainly the procedure of Gull 
& Daniel (1978) as described in Skilling & Bryan (1984). It determines the maximum 
entropy of an image consisting of model data points on a pre-set regular grid on condition 
that the model data, after convolving with the beam, fits the observed data according 
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to a chi-square fit. For large N, the program stops if x 2 is approximately equal to 
N + 3.29y/N where N is the number of observed data points (Skilling & Bryan 1984). 
This criterion corresponds to a 99 % confidence level for the chi-square fit. The observed 
data points were deconvolved and regridded with this program onto a 10" regular grid 
and smoothed with a 55" Gaussian beam. 

Taking the peculiar velocity structure of the Antennae (as seen, for example, in 
Ha (Amram et al. 1992)) into account and using spectral separation we can distinguish 
different components of the system that are not spatially resolved. Therefore we subdi- 
vided the total velocity range into two velocity ranges (v = 1420 - 1554 km s -1 and v = 
1550 - 1740 km s" 1 ) in order to distinguish between the interaction zone (at roughly v = 
1450 - 1590 km s" 1 ) and the nuclei (at v = 1600 km s" 1 and v = 1630 km s" 1 for NGC 
4039 (the southern galaxy) and NGC 4038 (the northern galaxy), respectively) (Amram 
et al. 1992, Rubin et al. 1970, Stanford et al. 1990). The cut in the velocity range is 
somewhat arbitrary since on the one hand the velocity gradient between the interaction 
zone and NGC 4039 is large, and on the other hand no clear spatial transition in the 
velocity occurs between NGC 4038 and the interaction zone. 

To obtain the [C II] integrated intensity maps in different velocity ranges we subdi- 
vided the raw data into two adjacent velocity bins and then applied the maximum entropy 
program for each velocity range separately. Since the upper and lower velocity ranges are 
derived only from part of the spectrum they have slighty different a values. The a values 
for the total, the upper, and the lower velocity ranges are 1.04 x 10~ 5 erg s~ 1 cm~ 2 sr~ 1 , 
0.77 x 10~ 5 erg s _1 cm _2 sr _1 , and 0.71 x 10~ 5 erg s _1 cm~ 2 sr _1 , respectively. 

4. Results 

4.1. Spatial Distribution of the [C II] Emission 

We superimposed maps of the total [C II] 158 /im line emission and the emission 
in the two seperate velocity ranges on an optical image (Laustsen, Madsen, & West 
1987) of the Antennae (Fig. la - c). For a better comparison we used the same contour 
levels in each velocity range. The contour levels are in steps of 1 a deduced from the 
data of the total velocity range. C + is seen over the full extent of the merging system 
(Fig. la). From the line intensity integrated over the total velocity range one can infer 
the origin of the bulk of the [C II] emission. In all of the three images the peak of the 
[C II] emission is situated at the interaction zone with slight differences in the location 
of the [C II] peak in each velocity range. The distribution of the [C II] emission in the 
total velocity range (Fig. la) and the upper velocity range (corresponding to the nuclei) 
(Fig. lb) are very similar and agree well with the overall shape of the galaxy pair as seen 
in short-exposed optical images, if we take the different beam sizes into account. The 
emission is elongated along the bridge between the two galaxies, including the interaction 
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zone, NGC 4039, and extends even further south. It is also extended in the north-west 
showing the western loop of [H II] regions around NGC 4038. However, for the upper 
velocity range one would have expected a stronger concentration of the [C II] emission 
at the nuclei and less strong at the interaction zone. Since the cut between the velocity 
ranges was made at 1550 km s -1 and the velocity range of the gas in the interaction 
zone is 1450-1590 km s" 1 (Amram et al. 1992; Rubin, Ford, & D'Odorico 1970; Stanford 
et al. 1990) it may be possible that most of the [C II] emission is coming from gas with 
velocities higher than 1550 km s _1 . This would also explain the relatively weak [C II] 
emission in the lower velocity range. The [C II] emission in the lower velocity range is less 
extended than in the other velocity ranges and more concentrated toward the interaction 
zone as expected from that velocity. The [C II] peak coincides also with the peak of the 
CO emission of the Antennae (Stanford et al. 1990) (Fig.2). Within the interaction zone 
the southern clump is the strongest CO emitting region. 



4.2. Estimate of the Parameters of the [C II] Emitting Region 

The results of the measurement of the [C II] 158 fim emission in the different velocity 
ranges for different positions are given in Table 1. The total [C II] luminosity from the 
region enclosed by the first contour line in Fig. la is L[cii] = 3.7 x 10 8 Lq for the total 
velocity range. In comparison with the [C II] ISO-LWS observations of the Antennae 
(Fischer et al. 1996), we determined a line flux of ~ 9 x 10~ 19 W cm -2 within the area 
of the larger ISO beam. This is about 2.4 times higher than the line flux measured 
with ISO. By weighting the data measured with FIFI within the ISO beam with a 
gaussion of 80" FWHM we obtain a line flux of 6.5 x 10 -19 W cm~ 2 which agrees within 
the statistical errors and the calibration errors with the value measured with ISO. A 
possible explanation for this discrepancy would be if the [C II] emission originates from 
small clumped regions spread over the interaction zone. Thus the emission of clumps at 
the edge of the ISO beam are only accounted for with the weight of the beam. 

The [C II] 158 /im fine structure line may arise from at least three components 
of the interstellar medium: photodissociation regions (PDRs), atomic gas clouds, and 
diffuse H II regions. In starburst galaxies and star forming regions the contribution of 
the [C II] fine structure line from PDRs is much stronger than from other components 
of the interstellar medium which can therefore be neglected. However, it is important to 
investigate the origin of the [C II] emission from the NGC 4038/39 system. We follow a 
procedure originally presented in Madden et al. (1993), to investigate the various origins 
of [C II] emission. 
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4-2.1. Estimate of the FIR continuum 

No FIR continuum map with a resolution comparable to the [C II] data exists for the 
Antennae. The only FIR continuum data available to compare with the [C II] emission 
is from IRAS observations. According to the formula given in Lonsdale et al. (1985) 

FIR = 1.26 x (2.58 x 1(T 14 x S m + 1 x 1(T 14 x S m ) [W m~ 2 ] 



we calculated the FIR luminosity from the IRAS flux densities at 60 /im and 100 fim 
for the entire Antennae system. With Sqo = 48.68 Jy and Sioo = 82.04 Jy (Surace 
et al. 1993), the FIR luminosity for the total system is Lfir = 3.6 x 10 10 Lq, with an 
assumed distance of 21 Mpc. Consequently the total [C II] luminosity is about 1% of 
the FIR luminosity. An attempt to obtain sufficient spatial resolution to resolve the 
nuclei from the IRAS data using HIRES and ADDSCAN was unsuccessful (Surace et al. 
1993). Therefore, to get an estimate of the FIR continuum separately for the interaction 
zone and the nuclei we use the FIR/radio correlation. However, in general the FIR/radio 
correlation is only valid for late-type galaxies. Despite this we use this correlation because 
the optical spectra of the nuclei indicate a mixture of H II regions with an old population 
and the interaction zone shows also a strong radio continuum (Hummel & van der Hulst 
1986) in combination with the young star forming region. The ISOCAM observations of 
the 6.7 and 15 /im continuum reveal very clumpy emission dominated by bright knots, 
and from the good correlation of the 15 /im continuum intensity with the [Ne III]/[Ne II] 
ratio, Vigroux et al. (1996) deduce that the 15 /im continuum is due to thermal emission 
from hot dust heated by the absorption of the ionizing photons emitted by young stars. 
On the kpc scale of our beam, the extension of the FIR continuum originating from 
thermal emission from moderately warm dust, and the 15 fim continuum created from 
these young star forming regions may be similar. Therefore, a strong but very confined, 
recent starburst as it is seen at the interaction zone with ISOCAM probably does not 
distort the FIR/radio correlation very much on our kpc scale. 

For the whole Antennae the FIR luminosity and the total flux density of the 20 cm 
radio continuum, Si.sghz = 486 ± 20 mJy (Hummel & van der Hulst 1986), fulfil the 
FIR/radio correlation: 

, -Pl.5GHz , £fir , , 

l0g WH^ = aXl0g ^T + 6 
with a = 1.30 ± 0.03, b = -25.82 and Pi. 5G h z = Si.sghz x 4nR 2 (Xu et al. 1994). 
Under the assumption that the FIR/radio correlation is also valid at the scale of our 
beam size we estimate the FIR continuum at the interaction zone and the nuclei by 
distributing the FIR continuum the same way as the low spatial resolution 20 cm radio 
continuum (Hummel & van der Hulst 1986) (Fig. ^). For that we convolved the 20 cm 
map to a spatial resolution of 55" and scaled it such that the total flux in the map 
matches the IRAS observation. We obtain a FIR continuum luminosity within our beam 
of L FIR = 2.0 x 1O 1O L at the [C II] peak, L FIR = 1.4 x 1O 1O L at NGC 4038, and 
L Fm = 8 x 1O 9 L at NGC 4039. 
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4-2.2. [C II] Emission From PDRs 

For the treatment of PDRs we used the model described by Stacey et al. (1991, and 
references therein). 

Photodissociation regions are the interfaces between H II regions and molecular 
clouds where photons with energies less than 13.6 eV escape from the H II regions, 
dissociate molecules and ionize elements with dissociation or ionization energies lower 
than the Lyman limit (13.6 eV). Carbon is the most abundant element with an ionization 
energy (11.3 eV) less than the 13.6 eV. Ionized carbon is then excited by collisions with 
electrons and/or atomic and molecular hydrogen. The gas in this region is heated mainly 
by photoelectric emission from grains illuminated by far-UV radiation. However, most 
of the far-UV flux which is absorbed by grains is converted into FIR radiation. 

Assuming high temperature (T ^> 91 K), high density (n ^> n cr it = 3.5 x 10 3 cm~ 3 ), 
optically thin [C II] emission, and a beam filling factor of unity we can estimate a lower 
limit of the C + column density and hydrogen mass (Crawford et al. 1985). Using a solar 
fractional carbon abundance of [C]/[H] ~ 3 x 10 -4 and assuming that all the carbon is 
in form of C + we derive a lower limit of the hydrogen column density. The result of 
this estimate is given in Table 2. For each of the positions we took the integrated [C II] 
intensity in the total velocity range. A comparison of the estimated masses of the PDRs 
with the molecular masses given by Stanford et al. (1990) show that the ratio of the mass 
of the PDR to the molecular mass is 6% for the whole Antennae, 6% for the interaction 
zone, 4% for NGC 4038, and 15% for NGC 4039. The high ratio for NGC 4039 is 
probably due to the fact that the beam of the [C II] observation also includes part of 
the interaction zone. Except for the high ratio in NGC 4039 the PDR to molecular gas 
mass ratio lies within the expected range of 1-10% (Stacey et al. 1991). 

To compare the [C II] integrated intensity with the CO (1 — > 0) intensity we use 
the CO-data from Aalto et al. (1995) since their observations were made with a single 
dish and with a similar beamsize, 43" FWHM. They observed three different positions: 
the interaction zone, NGC 4039, and NGC 4038. For the comparison, we use the [C II] 
emission in the total velocity range since the measurement of Aalto et al. (1995) was 
also obtained in a wide velocity range (~ 1300 km s _1 ). 

In Figure 4, we plotted Yjcii] = /[cii] /Xfir versus Y co = Iqo/Xfir- The advantage of 
taking these ratios is that they allow one to determine the gas density, n H , and UV-field, 
XuV) independently of the beam filling factor if it is assumed that the [C II], the CO, and 
the FIR emission arise from the same regions. The theoretically calculated ratios for a 
set of constant densities riu (solid lines) and a set of constant UV intensities xuv (6 eV 
< hv <13.6 eV) (thick dashed lines) are also shown in Figure ^ (Wolfire, Hollenbach , 
& Tielens (1989), Stacey et al. 1991). The observations for the interaction zone (filled 
triangle), for NGC 4039 (filled square) and for NGC 4038 (filled circle) are marked in the 
plot. The grey diagonal bars indicate the standard deviation of the FIR luminosity as 
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derived from the linear regression line of the FIR/radio correlation in Fig.l of Xu et al. 
(1994). We use this scatter as an estimate for the error in distributing the total IRAS 
flux over the mapped area according to the distribution of the 20 cm radio continuum 
map (see section 4.2.4). The results for the density and the UV-intensity are given in 
Table 3. Two possible solutions for riu and Xuv can be found: a low density (1) and a 
high density (h) solution. In both cases the UV-intensity is modest. Under the condition 
that xuv < 10 2 (low density solution) or Xuv/^h < 10~ 3 cm 3 (high density solution) a 
low /[cii]/^co ratio is predicted (Wolfire, Hollenbach, & Tielens 1989). In the low xuv 
regime the [C II] emission is reduced by low temperatures T < 92 K, and in the low 
Xvv/ n n regime by the CO self-shielding which brings the [C II] /CO transition closer 
to the surface and therefore lowers the C + column density. At the interaction zone a 
12 CO(2 - 1) / 12 CO(l - 0) line ratio of 1.2 has been observed and can be interpreted as 
arising from gas with n H > 10 4 cm -3 (Aalto et al. 1995, and Fig. 3 therein), supporting 
the high density solution in this region. For the nuclei this ratio has not been observed. 
The densities we derive for the high density solution are somewhat higher than what 
Wolfire, Tielens, & Hollenbach (1990) deduced for Orion, M82 and the inner region of 
the Galactic Center from the PDR model. However, we find a much lower FUV flux. 
The FUV flux for the high density solution is also at the lower limit of the range of the 
FUV flux derived from the ISO-LWS observations (Fischer et al. 1996). Fischer et al. 
(1996), however, determined a density for the PDRs which is an order of magnitude 
smaller than what we derived for the individual positions. 

We derived the beam filling factor for the single-component model (averaged over 
both components). To estimate the beam filling factor we used the ratio of the derived 
Xfir to the modelled xuv (Table 3). If we assume that the beam filling factor is unity 
and that most of the stellar photons are absorbed by dust grains, then xfir ~ 2 x xuv 
for B stars and xfir ~ Xuv for O stars. Including the beam filling factor, $, we get 
therefore xfir = (1 — > 2) x $ x xuv- For the individual galaxies NGC 4038 and NGC 
4039 where a large fraction of late type stars exists (Bushouse & Werner 1990, Keel et al. 
1985) the observed Xfir may not only be created by UV radiation, but also from visible 
light. Thus taking the beam filling factor as Xfir/xuv gives an upper limit. But there 
is still the uncertainty from the estimate of the FIR continuum. From the beam filling 
factors we estimate source sizes of ~ 24" for the interaction zone which is about the total 
size of the overlap region observed in CO (Stanford et al. 1990), and ps 19-21" for the 
nuclei which is about twice the size estimated from radio emission (Hummel & van der 
Hulst 1986) and from CO observations (Stanford et al. 1990). 

From the extinction corrected hydrogen recombination lines at the interaction zone 
observed with ISO (Kunze et al. 1996) we can also deduce the UV intensity originating 
in this region. An equivalent single star effective temperature of the stellar ionizing 
radiation field of 44000 K has been derived from the [Ne III] 15.6 /im / [Ne II] 12.8 fim 
line ratio (Kunze et al. 1996). This effective temperature corresponds to an 05 main 
sequence star. Therefore we use the properties of such a star to derive the Lyman 
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continuum emission. About 60 % of the total stellar luminosity is emitted in the Lyman 
continuum with an average photon energy of 18 eV (Panagia 1973). The rest is mainly 
emitted as far-UV radiation. Assuming case B recombination (Osterbrock 1989) we get 
for the relation between the extinction corrected Bra line intensity (Kunze et al. 1996) 
and the UV intensity 

Ajv = h-4 x x - — x [erg s cm z sr J 

/lZ/ 4 _ 2 J 5 -4 «4-2 

where is the relative intensity of Bra to H/3, a n-m is the total recombination 
coefficient, and 0:4-2 the effective recombination coefficient of H/3 given from Hummer 
& Storey (1987). From this expression we derive a UV intensity of I\jy = 345xo (for 
standard conditions: T = 10 4 K, n — 10 4 cm~ 3 ) and thus a far-UV intensity of Ifuv — 
230xo (xuv — 230) which is comparable to the far-UV we obtained from the ^jcii] — Yco 
plot for the interaction zone. We have done this estimate only for the interaction zone 
since ISO has only measured the hydrogen lines in the interaction zone. 



4.2.3. [C II] From Cold Neutral Medium (CNM) and/ or Warm Neutral Medium 

(WNM) 

In the neutral interstellar medium ionized carbon is excited by collisions with elec- 
trons and atomic hydrogen. The integrated intensity one would expect through collisions 
can be estimated by (see also Madden et al. 1993, 1997) 



hvA 
r[CII] = 17" 



f) 2exp(-f; 



1 + 2 exp(-f) + 1 + 2 exp(-f ) + ^ 



X c+ N u $ h 



where h is the Planck constant, v is the frequency of the transition, A is the Einstein 
coefficient for spontaneous emission, 2.29 x 10~ 6 s -1 (Nussbaumer & Storey 1981), T is 
the temperature, % is the density of atomic hydrogen, and X c is the ionization fraction 
of the medium. We assume the carbon abundance to be solar and that all the carbon is 
in the form of C + (Xq+ ~ 3 x 10~ 4 ). We used the peak value of the column density (Nn) 
plot of van der Hulst (1979, Fig. 5 therein): Ah ~ 6.3 x 10 20 cm~ 2 . The beam filling 
factor was assumed to be 1. The critical density for collisions with atomic hydrogen, 
^Hcrit, is deduced from the cooling function for collisional excitation of [C II] by atomic 
hydrogen (Launay & Roueff 1977), and that for collisions with electrons, n ecri t, is deduced 
from collision strengths for [C II] collisions with electrons (Blum & Pradhan 1992), and 
fitted as functions of temperature. 

The warm neutral medium in the Galaxy is characterized by low density, nu ~ 1 
cm -3 , and temperatures in the range T ~ 4 x 10 3 - 8 x 10 3 K, with an ionization fraction 
of X c ~ 3 x 10~ 2 (Kulkarni & Heiles 1987, 1988). The density in this medium is far below 
the critical densities for collisions of [C II] with electrons and atomic hydrogen. Thus 
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every collisional excitation will lead to radiative deexcitation. In this regime a maximum 
possible integrated [C II] intensity of I[cii] ~ 1 x 10~ 6 erg s _1 cm~ 2 sr~ 1 (at T = 4 x 10 3 K, 
nn = 1 cm -3 ) is estimated for collisions with both electrons and atomic hydrogen, 
and both collision partners contribute roughly equally to the integrated intensity. This 
estimated value is almost two orders of magnitude lower than the integrated intensity 
observed at the [C II] peak which coincides roughly with the peak of the H I column 
density. Since this estimated [C II] integrated intensity is also more than one order of 
magnitude smaller than the observed integrated intensity in the nuclei, where the H I 
column density should be smaller (van der Hulst 1979), a contribution of [C II] emission 
from the WNM to the observed integrated intensity is negligible. 

Carrying out the same calculations for cold neutral medium characterized by T ~ 50 
- 100 K, nu ~ 50 - 200 cm -3 , ionization fraction of X c m 5 x 10~ 4 (estimated for the 
Galaxy, Kulkarni & Heiles 1987, 1988) and again using the peak column density of 
A^h ~ 6.3 x 10 20 cm -2 (van der Hulst 1979) we get a maximum integrated intensity of 
J[cii] ~ 2.7 x 10 -5 erg s _1 cm~ 2 sr _1 (at T = 100 K, riu = 200 cm -3 ) and a minimum 
integrated intensity of I[cn] ~ 3.0 x 10~ 6 erg s" 1 cm _2 sr _1 (at T = 50 K, wh = 50 cm -3 ) 
for collisions with electrons and atomic hydrogen, where the contribution of collisions 
with electrons can almost be neglected. Therefore, at standard conditions (T = 70 K, 
n H = 100 cm~ 3 ) the [C II] emission from CNM is a factor of 9 smaller than the total 
[C II] radiation at the interaction zone, but may contribute as much as ~ 1 /s of the total 
[C II] emission. For the nuclei we estimate a similar contribution from CNM to the [C II] 
radiation using a H I column density of 3.5 x 10 20 cm -2 , from Fig. 5 of van der Hulst 
(1979) at the positions of the nuclei. 



4.2.4. [C II] From Ionized Gas 

Another possibility for the origin of the [C II] emission is from ionized gas, either 
extended low density warm ionized medium (ELDWIM) or "standard" H II regions. 

Although the shape of the contour plot of the [C II] fine structure line (of total 
velocity range) shows some similarities to the low resolution 20 cm radio continuum 
(Hummel & van der Hulst 1986) (Fig. 3), this extended radio continuum has a relatively 
steep radio spectrum (a < —0.8) hence it is almost completely non-thermal and does 
not trace extended ionized regions. The radio spectrum is steepest in the eastern part of 
the system, and regions devoid of strong optical emission, e.g. at the location of the dust 
patch in between the galaxy nuclei, show steeper spectra than regions that are bright in 
Ha and blue light. Therefore, the ELDWIM probably plays only a minor role and [C II] 
emission from this medium may be neglected. 

However, the radio knots 1-13 seen in the 4.9 GHz high resolution map of Hummel 
& van der Hulst (1986) (Fig. 5) have a sizeable thermal component. To estimate the 
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[C II] intensity expected to arise from the thermal emission of the knots within our beam 
we weight them relative to their distance from the center of the beam. We calculated the 
emission measure (Spitzer 1978) to be EM = 1250 cm _6 pc for the interaction zone and 
EM = 270 cm~ 6 pc and EM = 690 cm" 6 pc for NGC 4038 and NGC 4039, respectively. 
Making the crude approximation that all of the carbon in the H II region is in form of 
C + and that all of the observed EM results in the excitation of [C II], the expected [C II] 
integrated intensity is 



Here g u /gi = 2 is the ratio of the statistical weights in the upper and lower level and 
^crit = 49 cm -3 is the critical density of C + for collisions with electrons at 10 4 K (Blum 
& Pradhan 1992). Assuming an electron density of n e = 300 cm -3 as estimated from the 
extinction corrected ratio of the [S III] 18.71/zm/33.48/im fine structure line (Kunze et al. 
1996) the expected [C II] intensity is ~ 5.6 x 10~ 6 erg s _1 cm _2 sr _1 at the interaction 
zone, which is about a factor 20 smaller than the observed [C II] intensity in the lower 
velocity range. For the nuclei we assume an electron density of n c = 100 cm -3 (Rubin 
et al. 1970) and obtain expected [C II] intensities of I\cm = 3.3 x 10~ 6 erg s _1 cm _2 sr~ 1 
and I[cn] = 8.4 x 10~ 6 erg s^cnrtr- 1 for NGC 4038 and NGC 4039, respectively. This 
is about a factor 17 lower for NGC 4038 and a factor 8 lower for NGC 4039 than the 
observed [C II] intensities at the upper velocity range. However, most of the carbon in 
H II regions is probably in the form of C ++ rather than C + . Therefore, the contribution 
of the [C II] emission expected from H II regions is even lower and should be only a 
minor fraction of the observed [C II] radiation. 



As shown in the previous section a possible contribution of [C II] emission from the 
WNM, ELDWIM and H II regions to the total observed [C II] emission is small. For 
high temperatures and high densities a contribution of [C II] emission from CNM may 
become important. However, since at standard conditions of the CNM a contribution of 
[C II] emission to the observed [C II] emission is also small we will only consider PDRs 
as the origin of the [C II] emission in the following discussion. 

Stacey et al. (1991) demonstrated the utility of the [C II]/CO(1-0) ratio to dis- 
tinguish between starburst activity in galaxies and more quiesent regions. They find 
[C II]/CO ~ 6100 (4000)Q for starburst nuclei and Galactic starforming regions while the 
more quiesent regions show a ratio of ~ 2000 (1300) 1 . We find [C II]/CO ~ 2350 (1650)|] 





hvA 



Mil 

fli 



5. 



Discussion 



1 CO data (NRAO) corrected for main beam efficiency (~ 0.65) 
2 CO data (SEST) corrected for main beam efficiency (~ 0.7) 
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toward the [C II] peak (overlap region) and 1650 to 1800 (1150 to 1250) 2 toward the 
nuclei, implying there is no strong starburst activity taking place, given the resolution of 
our data. These ratios are more consistent with quiescent spiral galaxies, e.g. NGC 891 
and NGC 3628 (Stacey et al. 1991). However, observations with ISO (Kunze et al. 1996; 
Fischer et al. 1996; Vigroux et al. 1996), CO observations (Stanford et al. 1990; Aalto 
et al. 1995), and observations of the radio continuum (Hummel & van der Hulst 1986) 
indicate ongoing strong star formation activity in the interaction zone. A comparison of 
the CO flux of the interferometric observation (scaled to the size and the position of the 
single dish observation) of Stanford et al. (1990) (Fig. 2) with the single dish observation 
of Aalto et al. (1995) show a factor of about 4 more CO flux in the single dish observa- 
tion (Fqq) for the interaction zone and NGC 4039 and a factor 3 for NGC 4038. If we 
assume that l / 4 of the single dish CO flux at the interaction zone arises from a confined 
starburst region, F%% = Fgo x y 4 , and use a ratio of Ff^/Fgg = 6000 (the mean 
ratio observed for starburst galaxies, Stacey et al. 1991) for this region we find a ratio 
of (Fgfq - i^ciqV^co - F co) ~ 1200 for an underlying component. For F^ u] /F§ we 
used the observed ratio of 2350. This estimate suggests that the [C II] emission toward 
the interaction zone can also be explained by a two-component model with a confined 
starburst region and a quiescent surrounding molecular cloud system. Making the same 
estimate for NGC 4039 we get a [C II] to CO ratio of 400 for an underlying component. 
To make the same estimate for NGC 4038 one has to assume a lower [C II] to CO inte- 
grated intensity ratio for the confined starburst region since otherwise the ratio of the 
underlying component becomes negative. The lowest value for the [C II] to CO ratio in 
Stacey et al. (1991) is 340 for the SgrA +20 km s™ 1 cloud. If we assume a lower limit for 
the [C II] to CO ratio for the underlying component of 300 we find a maximum [C II] to 
CO ratio for a confined starburst region of 4400 for NGC 4038 using the observed ratio 
of FtQ^J Fqq = 1650. From this estimate we find that an enhanced star forming activity 
in NGC 4039 and a moderate activity in NGC 4038 could be explained by the data. This 
interpretation would be consistent with the results of the mid-infrared observations of 
NGC 4038/39 with ISOCAM (Vigroux et al. 1996). They show that the most active star 
formation in the Antennae system occurs in the overlap region, and within that region 
in a confined area (knot 2 in Fig. 5; knot A in Fig. 1 in Vigroux et al. 1996) coinciding 
with the southern clump of the CO emission in Stanford et al. (1990). The ISO SWS 
observations of the interaction zone can also be described by a young (ps 7 x 10 6 yr) 
starburst with an IMF extending up to 100 Mq (Kunze et al. 1996). The interpretation 
that the [C II] emission arises from confined small star forming regions is also consistent 
with the detection of Br7 knots in the interaction zone (Fischer et al. 1996). The case 
that most of the [C II] emission comes from confined star forming regions supports the 
PDR solution with a small beam filling factor and therefore with high density. The PDR 
solution with high density is also supported by the ISO observations of Fischer et al. 
(1996). 

Although we see moderate [C II] emission from the H II regions in the western loop 
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around NGC 4038, in the absence of CO data for that region and taking our beam size 
into account we cannot draw a conclusion about the star forming activity there. 

6. Conclusion 

We present a map of NGC 4038/39 in the [C II] 158 /im fine structure line. [C II] 
emission is detected over the optical extent of the system of galaxies and peaks at the 
interaction zone. The total luminosity of the [C II] line is 3.7 x 10 8 Lq which is about 
1% of the FIR luminosity of the Antennae. Only a negligible fraction of the observed 
[C II] emission can originate in the WNM, if conditions are similar to Galactic atomic 
clouds. Under normal conditions the [C II] emission from standard CNM makes only a 
small contribution to the total [C II] emission, however it may rise to Y 3 of the total 
[C II] emission for individual positions. Only minor fractions of the [C II] emission at 
the interaction zone and the nuclei can arise from H II regions. PDRs are the dominant 
source for the [C II] emission. We estimate minimum hydrogen masses associated with 
the [C II] emitting region of 1.9 x 10 8 Mq for the entire merging system and 6.8 x 10 7 Mq, 
3.2 x 10 7 Mq, and 3.7 x 10 7 Mq within one beam centered at the interaction zone, 
NGC 4038, and NGC 4039, respectively. Assuming a single emission component in 
the beam we derive a density of the [C II] emitting gas in PDRs of 1 x 10 5 cm -3 for 
the interaction zone and of 2 x 10 5 cm" 3 and 1 x 10 5 cm" 3 for NGC 4038 and NGC 
4039, respectively, and a far-UV intensity of 450x o for the interaction zone, 500x o for 
NGC 4038, and 250x o for NGC 4039. The derived beam filling factor of the emission 
from the PDRs is 20% for the interaction zone and 10-15% for the nuclei. However, 
the PDR model also allows a solution with low density (~ 10 2 cm" 3 ), high beam filling 
factor (~ 50%), and low FUV intensity (~ 100x o )- The low, beam averaged [C II] /CO 
ratio of 2350 toward the interaction zone and the even lower ratios at the nuclei indicate 
that no global starburst is going on either in the area surrounding the interaction zone or 
in the nuclear region. The high-excitation lines observed with ISO SWS which trace the 
starburst must therefore arise from a small, confined region in the interaction zone. This 
result is also supported from the observations with ISOCAM. Therefore on the scale of 
our [C II] beam a single emission component for the PDR is only a crude approximation. 
Using interfero metric and single dish CO observations and an expected [C II] /CO ratio 
for starburst regions and for quiescent clouds, we constructed a two-component model 
consisting of a confined starburst region and of molecular clouds enveloping the starburst. 
From this model we find that the [C II] emission at the interaction zone originates partly 
from confined starburst regions and partly from surrounding quiescent clouds. If we 
apply this model to the nuclei we get also enhanced star formation activity in NGC 4039 
with a low [C II] /CO ratio for the quiescent clouds but only moderate star forming 
activity in NGC 4038. This two-component model supports the high-density solution 
for the PDRs. 
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Future investigation of the Antennae in the [N II] 205 /xm fine structure line would 
be very helpful to further disentangle the origins of the [C II] line. Also observations at 
higher spatial resolution in the FIR regime (e.g. with FIRST and SOFIA) would be a 
great step forward to investigate this and other spatially very complex objects in more 
detail. 

We are grateful to the staff of the KAO for their competent support. This work 
was partially supported by the NASA grant NAG- 2-208 to the University of California, 
Berkeley. N.G. was supported in part by a Feodor-Lynen-fellowship of the Alexander 
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Table I. Integrated intensity in [CII] compared with CO (1 — > 0) data 



Source 


I[C1I] 
[10~ 5 erg s" 1 
cm^sr -1 ] 
(1420-1740 km s" 1 ) 


{[cii] a 

[10 5 erg s 1 
cm _2 sr _1 ] 
(1550-1740 km s" 1 ) 


/[Gil] a 
[10" B erg s" 1 
cm- 2 sr _1 ] 
(1420-1554 km s" 1 ) 


Ico b 

[10~ 8 erg s _1 
cm~ 2 sr _1 ] 


/[cii]/-fco c 


[C II] peak 


11.7 


6.8 


3.8 


4.9 d 


2350 


(Interaction zone) 








(7.0) e 


(1650) e 


NGC 4038 


5.5 


5.0 


1.3 


3.4 


1650 










(4.8) e 


(1150) e 


NGC 4039 


6.3 


3.9 


2.6 


3.5 


1800 










(5-0) e 


(1250) e 



a Notc that there is some overlap between the two subdivided velocity ranges due to the limited spectral resolution. 

b CO (1 — > 0) measurements of Aalto et al. (1995), scaled to the same units as [C II] intensity: 
1 K km s _1 = 1.6 X 10 -9 erg s -1 cm _2 sr _1 . 

C /[ CH ] of the total velocity range (1420-1740 km s _1 ) 

d the position of the beam of the interaction zone of Aalto et al. (1995) coincides with the [CII] peak in the total 
velocity range 

c CO data corrected for main beam efficiency (0.7) 
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Table 2. Minimum hydrogen mass and column density 



Source 


M min (H) 
1O 7 M 


-^min(Hj 

10 20 cm" 2 


[C II] peak 


6.8 


2.3 


(Interaction zone) 






NGC 4038 


3.2 


1.1 


NGC 4039 


3.7 


1.3 


total galaxy pair 


18.8 


0.67 (average) 
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Table 3. PDR model output 



Source 


PDR- 
solution a 


n H 
cm -3 


Xfuv b 


$b c 


Interaction zone 


1 


4 x 10 2 


150 


0.6 


([CII]) peak 












h 


1 x 10 5 


460 


0.2 


NGC 4038 


1 


1.4 x 10 2 


120 


0.50 




h 


2 x 10 5 


500 


0.1 


NGC 4039 


1 


3 x 10 2 


80 


0.45 




h 


1 x 10 5 


240 


0.15 



a l: low density solution 
h: high density solution 

b In units of Xo = 2 x 10~ 4 erg s" 1 cm" 2 sr™ 1 
(Draine 1978) 

c Beam filliner factor: $ h = 2mK - 
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Fig. 1. — Integrated [C II] 158 /iin intensity map in the velocity ranges of a) 1420- 
1740 km s _1 , the total [C II] emission, b) 1550-1740 km s™ 1 , the velocity range contain- 
ing the 2 nuclei, and c) 1420-1554 km s _1 , the velocity range of the interaction zone, 
superimposed on an optical image of NGC 4038/39 (Laustsen, Madsen, & West 1987). 
The contour levels in all three maps are the same, starting at 1 a with steps of 1 a 
(1.04 x 10" 5 erg s~ 1 cm~ 2 sr~ 1 ). The hatched circle indicates the beam (FWHM of 55"). 
The symbol "S" indicates the southern clump of the interaction zone in Stanford et al. 
(1990). 



Fig. 2. — CO (1-0) map of Stanford et al. (1990) (white contours) superimposed on the 
[C II] integrated intensity map of the total velocity range (black contours). The FWHM 
of the [C II] beam, indicated by the hatched circle, is 55". S, W, N, E, indicate the 
southern, western, northern, and eastern CO clump at the interaction zone (Stanford 
et al. 1990). The grey area represents the area covered by the [C II] arrays. (Note: The 
levels of the CO contours at NGC 4038 are different from the levels at the interaction 
zone and NGC 4039) 



Fig. 3. — 1.5 GHz continuum map of Hummel & van der Hulst (1986) (white contours) 
superimposed on the [C II] integrated intensity map of the total velocity range (black 
contours). The FWHM of the [C II] beam, indicated by the hatched circle, is 55"; the 
HPBW of the radio map is 16". The grey area represents the area covered by the [C II] 
arrays. 



Fig. 4. — Ratios ^jcii] — -^[ciii/xfir vs. Yco — -^co/xfir- The solid lines represent a set 
of constant densities and the dashed lines a set of constant UV intensities characteristic 
of PDRs (Stacey 1991). The filled triangle marks the observations at the [C II] peak, the 
filled square NGC 4039, and the filled circle NGC 4038. The grey diagonal bars indicate 
the standard deviation of xfir- 



Fig. 5. — 4.9 GHz continuum map of Hummel & van der Hulst (1986) (white contours) 
superimposed on the [C II] integrated intensity map of the total velocity range (black 
contours). The FWHM of the [C II] beam, indicated by the hatched circle, is 55"; the 
HPBW of the radio map is 6". The numbers 1-13 indicate the discrete radio knots 
(Hummel & van der Hulst 1986). The grey area represents the area covered by the 
arrays. 
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